3-point temperature anisotropies in WMAP: 
Limits on CMB non-Gaussianities and non-linearities 
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We present a study of the 3-pt angular correlation function u>3 =< 816283 > of (adimensional) tem- 
perature anisotropies measured by the Wilkinson Microwave Anisotropy Probe (WMAP). Results 
can be normalized to the 2-point function u>2 =< 8182 > in terms of the hierarchical: q-j, ~ W3/W2 or 
dimensionless: dz ~ W2/wfJ 2 amplitudes. Strongly non-Gaussian models are generically expected 
to show ds > 1 or q$ > 10 3 d3. Unfortunately, this is comparable to the cosmic variance on large 
angular scales. For Gaussian primordial models, 53 gives a direct measure of the non- linear correc- 
tions to temperature anisotropies in the sky: 8 = 8l + /jvlt(<5|,— < 8\ >) with /jvlt = 93/2 for 
the leading order term in wf. We find good agreement with the Gaussian hypothesis ds ~ within 
the cosmic variance of ACDM simulations (with or without a low quadrupole). The strongest con- 
straints on §3 come from scales smaller than 1 degree. We find §3 = 19 ± 141 for (pseudo) collapsed 
configurations and an average of 93 = 336 ±218 for non-collapsed triangles. The corresponding non- 
linear coupling parameter, /jvl, for curvature perturbations $, in the Sachs- Wolfe (SW) regime is: 
fxL = 93/6, while on degree scales, the extra power in acoustic oscillations produces /jvl ~ 93/30 
in the ACDM model. Errors are dominated by cosmic variance, but for the first time they begin to 
be small enough to constrain the leading order non-linear effects with coupling of order unity. 
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I. INTRODUCTION 

We study the cosmic microwave background (CMB) 
anisotropies measured by WMAP |l| from the point of 
view of the 3-point angular correlation W3 . Current mod- 
els for structure formation predict that the initial seeds, 
that later give rise to the observed structure in the uni- 
verse , were Gaussian. If so, one would expect the angu- 
lar 3-point CMB temperature correlations to be almost 
zero. Even for a Gaussian field, any non-linear process- 
ing of the primordial perturbations would lead to small 
non-Gaussianities. For an excellent review see and 
references therein. 

Angular correlations have some advantages over spher- 
ical harmonic analyses. The main disadvantage is that 
different angular bins are highly correlated, which means 
that we need to use the full covariance matrix to assess 
the significance of any departures between measurements 
and models. The advantage of using the N-point function 
is that correlations can be easily sampled from any re- 
gion of the sky, even when a foreground mask is required, 
provided the region is large enough. In comparison, har- 
monic decomposition breaks the angular symmetry in the 
sky into orthogonal bases which have an arbitrary phase 
orientation. This is not a problem in the case of full sky 
coverage, but the results will be coordinate-dependent if 
a mask is used. Moreover, higher order correlations, such 
as the bispectrum, have to focus on some small portion of 
the mutidimensional space of triangle configurations (see 
0)- The formal equivalence between the bispectrum and 
the 3-point function is hardly realized in practice and as 
such it is of interest to study both separately. 

Our analysis procedure is similar to that presented 
in |3J|, to which the reader can refer for further de- 



tails. We use the same simulations where we've gener- 
ated maps using the standard ACDM model and a low-Q 
ACDM model, a fiducial variation of the ACDM model 
with low values of the quadrupole and octopole. We 
concentrate our analysis on the V-band WMAP maps 
(FWHM = 2T at 61 GHz) generated using HEALPix 
with pixel resolution n=512 (7' pixels) and to which we 
apply the conservative WMAP kpO foreground mask 0] , 
The resulting maps have approximately 2.4 million pix- 
els. 



II. CORRELATION FUNCTIONS FROM WMAP 

The 2-point angular correlation function is defined as 
the expectation value or mean cross-correlation of tem- 
perature fluctuations: 
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at two positions qi and q2 in the sky: 
w 2 (9) = (<J(q 1 )«y(q 3 )) ) 
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where 9 = |q2 — qi|, assuming that the distribution is 
statistically isotropic. We normalize fluctuations to be 
dimensionless, To = 2.73K, or dimensional with To = 1. 
To estimate u>2{0) from the pixel maps we use: 
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where Si and Sj are the observed temperature fluctua- 
tions on the sky and the sum extends to all pairs i,j 
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separated by ± Ad. The mean temperature fluctuation 
is subtracted so that (Si) = 0. The weights Wi can be 
used to minimize the variance when the pixel noise is 
not uniform, however this introduces larger cosmic vari- 
ance. Here we follow the WMAP team and use uniform 
weights (i.e. Wi = 1). At the largest scales we use bins 
A0 whose size is proportional to the square root of the 
angle and use linear bins of A9 ~ 0.26 degrees at the 
smallest scales. 



A. 3-point function 

In a similar fashion, the 3-point angular correlation 
function is defined as the cross-correlation of temperature 
fluctuations at 3 positions qi, q2 and q3 on the sky: 



W3(9 12 ,9 13 ,9 2 3) = (S(qi)S(q 2 )S(q 3 )), 
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where 9ij = |qj — qj|. Consider a triplet of pixels with 
labels 1,2,3 on the sky. Let 6*12 and #13 be the angu- 
lar separations between the corresponding pairs of pixels 
and a the interior angle between these two sides of the 
triangle. Assuming that the distribution is statistically 
isotropic, W3 only depends on the sides of the triangles 
that define the 3 positions on the sky. One can char- 
acterize the configuration-dependence of the three-point 
function by studying the behavior of W3 (a) for fixed #12 
and #13. 

To estimate ^3(^12, 0i3, #23) from the pixel maps we 
use: 

W 3 (012, 6»i3, 023 = ^ ) (5) 

T,i,j,k "'< "'./ " 

where Si, Sj and 5k are the temperature differences in 
the map and the sum extends to all triplets separated by 
0i2 ± A0, 0i3 ± A0 and 023 ± A0. We use same bins and 
weights as in the 2-point function above. 

The 3-point function can either be normalized to a 
dimcnsionlcss amplitude: 
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where 



•Jij = w 2 (9ij), or a hierarchical amplitude: 
("3 



?3 



W 12 W 2 3 + W 12 Wi3 + W23W13 



Errors 



(6) 



(7) 



The covariance matrix between different angular bins 
can be calculated from the simulations: 



Ca = (Aw 2 {9 l ) Aw 2 (9 3 )) 
1 N 

= ^E AW 2(^)A W 2 L (^), 
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where Aw\ (0j) = w%{6i) - wi{9 t ). Here w% (0,) is the 2- 
point or 3-point function measured in the L-th realization 
(L = 1 . . . N) and 102(0,) is the mean value for the N 
realizations. The case i = j gives the error variance: 

An estimate of the errors for correlations measured on 
the real sky can be obtained using a variation of the 
jackknife error scheme. This has the potential advan- 
tage of producing an error estimate which is model inde- 
pendent. The accuracy of the jack-knife covariance have 
been tested for both WMAP Q and the APM and SDSS 
survey EJ, 0- In the estimation, the sample is first 
divided into M separate regions on the sky, each of equal 
area. The analysis is then performed M times, on each 
occasion removing a different region. These are called 
the (jackknife) subsamples, which we label k = 1 . . . M. 
The estimated statistical covariance for w 2 at scales 0, 
and 0j is then given by: 
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where 



i) is the measure 2 or 3-point function in the 



k-th subsample (k = 1 . . . M) and w 2 (9i) is the mean 
value for the M subsamples. The case i = j gives the 
error variance. Note how, if we increase the number of 
regions M, the jackknife subsamples are larger and each 
term in the sum is smaller. We typically take M = 8 
corresponding to a division of the sphere into 8 octants. 
Note that the 3-point correlation function is in general a 
function of 3 variables. We will fix two of these variables 
and show variations of w 3 or q 3 as a function of a single 
variable, for which we estimate the covariance with the 
above prescription. 



C. The algorithm 

Estimation of correlations in configuration space is 
quite time consuming. A brute force algorithm to es- 
timate all pair separations in a map with N pixels takes 
~ N 2 operations. To do all triplets in Eq.0 we need 
~ iV 3 operations. In the healpix n — 512 WMAP maps 
(7 arcmin pixel's) we have N ~ 2.4 x 10 6 (with the 
KpO mask galactic cut) so a brute force analysis requires 
~ 10 20 operations, well beyond current computing power. 
We overcome this difficulty in two steps. For large scales, 
> 1 degree, we reduce the resolution of the maps to 
n = 64 (1 degree resolution) by averaging fluctuations 
over nearby pixels. The resulting correlation function 
contains the same information at scales larger than the 
new resolution scale. Now the number of pixels is re- 
duced to N = 3 x 10 4 and the number of operations to 
~ 10 13 , which lies within current computer power (Ter- 
aflops). For scales smaller than a few degrees we can 
reduce the number of operations by precalculating the 
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FIG. 1: 3-point function W3 =< S1S283 > at large angular 
scales. Bottom spanel shows W3 for equilateral triangles as a 
function of side 612. The top panel shows results for isosceles 
triangles #13 = O12 with #23 = l.ldeg (long dashed line), 
823 = 812/2 (short dashed line) and 823 = 28\2 (continuous 
line), in all cases as a function of side 9\2- Shaded regions are 
for the 68% confidence levels estimated using jackknife errors. 



FIG. 2: Reduced 3-point function 53 at small angular scales. 
Different panels show configurations with triangles of fixed 
sides 812 and #13 (as labeled in each panel) and variable 823 
which is given by the interior angle a between #12 and 813 
(a — corresponds to 823 = and a — 180 corresponds 
to 823 = 201a). Shaded regions are for the 68% confidence 
levels estimated using low-Q ACDM simulations. Errorbars 
are from jackknife variance. 



distance between neighboring pixels. The number of op- 
erations is then reduced toJVx M 2 , where M << N is 
the number of neighbors. In our case we have M ~ 10 3 
for pixels around 9 < 4 degrees, and thus N x M 2 ~ 10 13 
as in the case for the larger angles. In our implementa- 
tion we have an overlap between the large and small scale 
estimations at around 2 — 4 degrees. The agreement be- 
tween both estimates here is quite good for both W2 and 
W3. In the case of the 2-point function we have compared 
both estimates in the whole range 2 — 180 degrees and 
find excellent agreement. 



D. Results 

Results are presented in Figures ^ and [21 At large 
scales Fig^ shows the results of W3 in units of (50fiK) 3 
(eg T = 50/i-K in Eq.P). Note that vj}J 2 ~ 50fiK at 
small scales, so that W3 ~ = ws/w^ 2 in these units 
(see Eq.jS]). Errorbars in both cases are from jackknife 
subsamples and agree well with errors estimated from 
the simulations (see Gaztanaga etal 2003 for more de- 
tails). The results for the equilateral and collapsed tri- 
angle configurations can be compared to the COBE anal- 
ysispresented in Fig. 1 of ; 8j . The pseudo-collapsed case 
in y| would correspond here to #23 = l.ldeg here. We 
can see how the COBE and WMAP results are in good 
agreement where the two curves even see the very same 



bumps and valleys. This is not totally surprising given 
the very good match between COBE and WMAP tem- 
peratures maps, pixel by pixel, when smoothed on large 
scales. The top panel in Figure ^ a lso shows two new tri- 
angle configurations: one that is elongated, 823 = 6*12/2 
(short dashed line), and a wide triangle, 6*23 = 2#i2 (con- 
tinuous line). The results for intermediate angular con- 
figurations are similar. In all cases the agreement with 
the Gaussian prediction W3 = is excellent, given the 
strong covariance. 

At small scales we show the 3-point function in terms 
of the reduced amplitude 173 in Eq. [7] where the best con- 
straints are found on the smallest scales. As can be 
seen in the figure, as we approach triangles with side 
6 ~ ldegree, the cosmic errors become large Aq3 ~ 1000. 
We have fit the data to a constant value of 93, indepen- 
dent of scale using a % 2 -test: 



(ii) 



where Aj = q^icti) — q§ is the difference between the 
"observation" O and the fitted value. In order to elim- 
inate the degeneracies in the covariance matrix, we per- 
form a Singular Value Decomposition (SVD) of the co- 
variance matrix (see Gaztanaga et al. 2003 for more de- 
tails). A joint fit with the jackknife covariance matrix for 
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collapsed configurations with 9 < 1 degree, yields: 

q coUapsed = jg ± m 

The best signal-to-noise is found for the non-collapsed 
triangle configurations with 9 ~ 0.5 degree we find: 

q 3 = 336 ±218 ; 9 = 0.5 cfep (13) 

providing a marginal detection at about 90% confidence 
(shown as dashed line in FigJ2J. This fit corresponds 
to the jackknife covariance matrix. A similar result is 
found for the covariance matrix in the low-q ACDM sim- 
ulations. The errors increase by almost a factor of two 
when using the covariance of the standard ACDM sim- 
ulations. This increase is due to the larger cosmic vari- 
ance introduced by the quadrupole and octopole which 
are larger than in the WMAP data (see 0, citepGaz03 
for a discussion). 

III. DISCUSSION 

The above results for the 3-pt function confirms previ- 
ous analyses (Q and references therein) that find WMAP 
to be in good agreement with the Gaussian hypothesis. 
Strong non-Gaussian statistics with d 3 in Eq.jS] of order 
unity or larger are ruled out. On larger (COBE) angular 
scales the results are dominated by cosmic variance so 
that A<i 3 ~ 1 0| and Aq 3 ~ 1000. On subdegree scales 
errors are smaller and we find evidence for a marginal de- 
tection (90% level) of non-Gaussianity in Eq.^31- This 
result can be taken as a detection or, at a higher sig- 
nificance level, a bound on possible non-linear effects. 
Expanding the observed temperature fluctuation S as a 
perturbation over some primordial Gaussian linear con- 
tribution 5^: 

S = S L + f NLT (5l-<S 2 L >), (14) 
we then have that to leading order in W2- 

W 3 = 2 f NLT (W12W23+W 12 W1 3 +W23W 13 ) + 0[W2} (15) 



so that q 3 ~ 2 Jnlt- The corresponding non- linear cou- 
pling parameter, /nl, for curvature perturbations $ is: 



* = $l + /m ($!-<$! >). (16) 

In the Sachs- Wolfe (SW) regime (S ~ <I>/3) we have: 
fxY — 93/61 while on degree scales, the extra power 
in acoustic oscillations produce //vl — 93/30 in the case 
of ACDM 11]. This corresponds to f NL ~ 11.2 ± 7.3 
for Eq.^21 and /jvl — 0.6 ± 4.7 for the collapsed case 
Ea. 1121. For a comparison with recent predictions see 
UllJlUQEl and references therein. 

Our bounds on Jnl seem marginally better than the 
ones presented in [2J from the WMAP bispectrum, but a 
direct comparison is not straight forward. First of all, our 
modelling provides a measurement of non-linearities in 
the final AT while Q provide a modelling for non-linear 
primordial curvature perturbations before radiation 
transfer. Thus while limits on q 3 are model independent, 
direct limits on fpj l depend on the assumed model for ra- 
diation transfer. Secondly, the bispectrum configurations 
used by Q are a subset of all possible triangle configu- 
rations so that both approaches arc neither equivalent 
nor optimized in the same way. In both cases, errors 
are dominated by cosmic variance, and they begin to be 
small enough to constrain the leading order non-linear 
effects with coupling /jvl which are of order unity. 
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